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y—i .We clarify some misunderstandings on the time-dependent current density functional theory for open quantum systems we have 
O "recently introduced [M. Di Ventra and R. D'Agosta, Phys. Rev. Lett. 98, 226403 (2007)]. We also show that some of the recent 
formulations attempting to improve on this theory suffer from some inconsistencies, especially in establishing the mapping 

Q_i between the external potential and the quantities of interest. We offer a general argument about this mapping, showing that it 

& must fulfill certain "dimensionality" requirements. 
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Experience teaches us that no physical system can be con- 
sidered as completely closed, in the sense that some degrees 
of freedom are either too many to consider explicitly, or their 
microscopic features are not of interest, and only their macro- 
scopic thermodynamic properties are of importance. These de- 
grees of freedom are then treated as baths or reservoirs of the 
system of interest, making the latter "open". Nonetheless, the 
theory of open quantum system is relatively new, starting with 
the pioneering works by Zwanzig and Nakajima, 12 who worked 
out the theory in terms of the statistical operator (density ma- 
trix), and more recently with the formulation of the theory of 
open quantum system in terms of a state vector ("wavefunc- 
tion")..M 

In parallel, another theory, Density Functional Theory (DFT), 
developed since the first papers by Kohn and co-workers, 5,6 has 
changed the way we think in terms of many-body quantum me- 
chanical systems. Indeed, the successes of DFT are countless, 
and the theory has now become a standard for many calcula- 
tions of electronic structure, electron transport, materials prop- 
erties, and chemical reactions.— ~— 

The attempt to bring together the two theories appears there- 
fore completely natural. 13-19 In particular, the present authors 
have proved a theorem that establishes a one-to-one correspon- 
dence between a vector potential acting on an open quantum 
system and its current density a theorem that was later ex- 
tended to the correlated motion between electrons and (in prin- 
ciple quantum) ions, 16 2() This theory, which we named stochas- 
tic time-dependent current-DFT (STDCDFT), then allows for 
the solution of a many-body open quantum system with effec- 
tive single-particle equations, a significant computational sim- 



a ETSF Scientific Development Center, Departamento de Fisica de Materiales, 

Universidad del Pais Vasco, E-20018 San Sebastian, Spain 

b IKERBASQUE, Basque Foundation for Science, E-48011, Bilbao, Spain 

c Department of Physics, University of California - San Diego, 92093 La Jolla, 

USA 

\ Email: roberto.dagosta@ehu.es, \ Email: diventra@physics.ucsd.edu 



plification. 

Although still in its infancy, STDCDFT for open quantum 
systems has received a certain amount of attention, especially 
for its promise to expand the field to systems previously in- 
accessible within standard DFT. However, there is also some 
confusion regarding its foundations which is mostly evident in 
recent developments, for instance a new result that shows how 
in principle we can mimic the dynamics of an open quantum 
system by reverting to the analysis of a closed quantum system 
made of non-interacting particles. 1718 This result would fur- 
ther dramatically reduce the computational cost for evaluating 
the dynamics of the open quantum system. 

The reason for this reduction is simple. In general, the 
open quantum system could be described in two ways. On the 
one hand, we can use the many-body statistical operator, or the 
single-particle reduced density matrix. However, for large sys- 
tems, the number of degrees of freedom, i.e., the dimensional- 
ity of the operator we need to consider, scales badly although 
some schemes can be devised to cure this scaling. 21 On the 
other hand, a formulation in terms of a stochastic state vec- 
tor is computationally cheaper than a density matrix one, but 
nonetheless it requires the consideration of a statistical ensem- 
ble of replicas of the system. Larger the ensemble, finer will 
be the physical information we can extract from the dynamics. 
Instead, if we could describe the dynamics of the open system 
with a closed non-interacting system, we would have reduced 
a complex and expensive problem to a size similar to that we 
usually tackle with our present computational means. 

Obviously,-a situation common in DFT-we cannot expect 
to have access to the complete physical information about the 
open quantum system. We must accept that our non-interacting 
doppelganger will be able to describe only some quantity of 
choice, normally either the single-particle density or the single- 
particle current density. In principle, any other physical ob- 
servable can be extracted provided we know its expression in 
terms of these "fundamental" quantities. However, in practice, 
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we rarely know how to write such obervables in terms of the 
density or the current density, thus limiting the range of physi- 
cal quantities that can be actually computed. Our goal for this 
paper is to both clarify some of the theoretical foundations of 
STCDFT that do not seem to have been appreciated in the lit- 
erature, as well as to show that some of the recent attempts to 
extend this theoryi^l^ suffer from inconsistencies, in particular 
those related to the mapping between external potentials and 
the basic quantities of any density functional theory. 

Let us then start by noting that the choice about which quan- 
tity we want to base our theory on is usually made a priori. 
For this reason we usually talk about DFT or Current-DFT.^S 
For closed systems, Current-DFT is regarded as a more com- 
plete theory than DFT since, via the continuity equation, one 
can calculate the single particle density starting from the single 
particle-current density. The situation is instead more complex 
in open quantum systems since, as we are going to show in 
the following, prima facie the continuity equation would seem 
not to uniquely determine the single-particle density, given the 
single-particle current density. 

Let us then begin by recalling that in the theory of open 
quantum systems, the equation of motion for the ensemble- 
averaged particle density, n(r, t), is given by (see, e.g.,—) 



d t n(r,t) = -V-j{r,t) + T B (r,t) 



(1) 



where Ts{r : t) describes the density modulation induced by 
the presence of the bath(s) and j (r, t) is the ensemble-averaged 
current densitjQ. Eqn[T]can be obtained from the master equa- 
tion of motion for the density matrix, p (ft = 1) 



dtp(t) 



H(t),p(t) 



dt'K(t,t')p(t'). (2) 
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Here, K (t, t') is a memory kernel that describes the action of 
the bath(s) on the system-provided it is chosen to preserve the 
positivity of the density matrix at any given time 'iW 4 - 25 and 
H(t)is the Hamiltonian of the system that evolves under the ac- 
tion of an external vector potential (possibly time dependent), 
A(r, t), and in the presence of a particle-particle interaction op- 
erator W, 
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We work here in a gauge where the scalar potential has been set 
to zero at all times. 

In general, little is known about the term Tb appearing 
in the continuity equation ([T). Although not strictly neces- 
sary for the considerations we make later, it would be useful 
if Tb could be written as the divergence of some current. This 

* In this theory we do not allow for the system to exchange particles with the 
environment. 



would immediately imply the conservation of the total num- 
ber of particles. Recently, Gebauer and Car have shown that 
this is indeed the case in the Markov approximation for cer- 
tain couplings between the system and the environment ,2&£L 
More recently a formal condition on current conservation and 
its use to test the validity of the approximations made to ar- 
rive at the master equation has been proposed.— If Tb is the 
divergence of some current-this extra current being generated 
by the exchange of momentum and energy with the external 
environment-then the continuity equation assumes the usual 
physical meaning and particle number conservation is guaran- 
teed. Moreover, one could assume that 



T B = -V • 



eC(r, t) 



n(r,t) 



(4) 



where without loss of generality we have written the extra cur- 
rent as en(r, t)C(r, t) / m. If this were the case, the total current 
3t(?i t) = j{ r i t) + en(r, t)C(r, t)/m fulfills the standard con- 
tinuity equation, 



d t n(r,t) = -V-j T (r,t). 



(5) 



In the total current, C (r, t) appears then to play the role of an 
extra vector potential, effectively inducing a diamagnetic cur- 
rent in the system, en(r,t)C(r,t)/m. Interestingly, in Ref.— 
C(r, t) has been dubbed "leakage potential", although we would 
like to point out that its definition from eqn |4] essentially pre- 
vents any leakage, i.e. loss of particles, from the system. 

In standard TDCDFT for closed quantum systems, where 
the extra term in Eq. (Q3 is not present, i.e., Ts{r : t) = 0, 
the continuity equation establishes an important link between 
the single -particle current density j(r, t) and the single -particle 
density n(r, t). Indeed, this equation allows for expressing the 
latter in terms of the former via a simple time integration. This 
seems a trivial point, but let us clearly state it: The continu- 
ity equation, when T B (r,t) — can be used together with 
the initial conditions to uniquely determine the single-particle 
density given the single-particle current density. The basic the- 
orem of TDCDFT therefore establishes a one-to-one mapping 
between the current density and the vector potential applied to 
the system, once the initial conditions are provided. 10 i 12 i 29 In 
fact, whenever one discusses the uniqueness of the mapping 
between the single-particle density and single-particle current 
density with the external vector potential, it is implicitly un- 
derstood that the particle density is obtained from the current 
density via the continuity equation. 

The situation becomes a bit trickier whenever Tb (r, t) ^ 0. 
In this case one needs the expressions of both j and Ts{r, t) to 
be able to determine the single-particle density. However, writ- 
ten in this way the problem corresponds to the solution of a lin- 
ear differential equation with non-constant coefficients. Then, 
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if one knows Tb we just need to perform a simple time integra- 
tion. From a TDCDFT formulation of the theory then, as it has 
been shown,i^iI£ one can build a one-to-one mapping between 
the (average) current density and the external vector potential, 
provided the initial condition and the bath operators are fixed. 

Here we stress that the inclusion in the mapping of the den- 
sity along with the current density-as for example it has been 
tried in Ref.i^-is not just redundant, but simply incorrect. In- 
deed, we will show in the following that the mapping between 
the set of current and particle densities and the vector potential 
is not unique. Our argument can be extended to other cases. 
For example, along the same lines we can prove that the map- 
ping between the current density and the scalar potential can- 
not be one-to-one. One of us (RD) has already presented an- 
other proof of the same result, in particular showing that the 
current density is in general non-V-representable. 30 This, there- 
fore, precludes the access to all components of the current den- 
sity by standard TDDFT. Our following argument also sets a 
necessary condition for building any density functional theory, 
in terms of establishing a one-to-one correspondence between 
two physical quantities. This condition states that the two quan- 
tities must have the same "dimensions": a vector with a vector, 
a scalar with a scalar and so on. Any other mapping will be 
either redundant or wrong. 

That this point is still not clear in the community is evident 
for example in the thesis of Theorem 1 of Ref. 17 where it is 
stated that there exists a one-to-one correspondence between a 
vector potential and the set of functions given by the current 
and particle densities. The Authors of that work go also as far 
as clearly stating in their discussion of the content of Theo- 
rem 1 that their aim is to consider the particle density on the 
same footing as the current density. However, this statement 
cannot be correct as a simple mathematical characterization of 
the functional spaces connected by this mapping shows. If the 
particle and current densities were independent functions, then 
they would belong to a 4-dimensional functional space. This 
amounts to saying that Tb (r, i) is not completely determined 
by the sole knowledge of n(r, t) and j(r, t). Indeed, if Tb{t, t) 
was determined by the sole knowledge of the single-particle 
densities, then we could use the continuity equation (fl~|i as a 
non-linear differential equation for the density and by solving 
it, obtain n(r,t) in terms of the initial conditions and j(r,t). 
Since a general vector potential is a vector of three independent 
functions, Theorem 1 in 17 would imply that a 4-dimensional 
functional space formed by the density plus the three compo- 
nents of the current density is locally homeomorphic to a 3- 
dimensional functional space spanned by the vector potential. 
This is obviously incorrect. 

The solution to this inconsistency is that the continuity equa- 
tion (Q~|) is indeed a non-linear equation in the particle and cur- 
rent densities from which the particle density can be obtained 
starting from the current density, once the bath operator and ini- 



tial conditions have been fixed.^^ This amounts to saying that 
J 7 b(i'> t) is a functional of n(r, t) and j(r, t), or better of j(r, t) 
alone, and that equation ((TJ admits a unique physical solution. 

In fact, we can provide a simple counter argument to the ex- 
istence of a mapping between the vector potential and the set of 
functions given by the particle and current densities. Consider 
the extremely simple case in which the response functions of 
the system are local in time and constant in space, and that the 
(full) response is linear. This case is obviously extremely sim- 
plified. However, to provide a counterexample to a theorem we 
do not need to discuss more complex cases. For this oversim- 
plified system we can write the following response equations 
(for simplicity, in the next few equations we drop the spatial 
and temporal dependence of j and n) 



( jx 

jy 

jz 

V n J 



= r 




(6) 



where T is a 4 x 3 matrix of constant coefficients of the system 
response. The existence of a one-to-one mapping between the 
particle and current densities and the vector potential would 
then suggest that we can write a similar equation, starting from 
the particle and current densities 
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(7) 



where now V is a 3 x 4 matrix. Combining these results it 
follows that 

ff ' = U; f'f = 1 3) (8) 

where 1^ is the N x N identity matrix. 

We can now easily prove that either the current and parti- 
cle densities are not independent (and therefore the mapping is 
indeed redundant), or the mapping between the current and par- 
ticle densities and the vector potential is not one-to-one. Due to 
the particular structure of eqn|6] we can invert the relations be- 
tween the current density and the vector potential and express 
the latter in terms of the former. Since this is in principle a 3- 
by-3 system of linear equations, it admits always at least one 
solution. Then, we can use the results found in this way in the 
equation for the particle density, therefore expressing the latter 
in terms solely of the current density, thus proving the mapping 
redundant. 

On the other hand, since eqn (0 is a system of 3 equa- 
tions in 4 variables, we can always find two distinct sets of 
current and particle densities providing the same vector poten- 
tial. Therefore an infinite number of solutions can be found. 
To show this we proceed in this way: Assume that jo, no are 



3 



a set of current and particle density, and that Aq is the result- 
ing vector potential obtained from (0. Now let us consider the 
particle density n = no + An. We want to show that we can 
find another current density j = jo + Aj which together with n 
gives, via (0, the same vector potential. Since An is arbitrary 
we can choose as initial condition An(r, t = 0) = 0, therefore 
the two sets of current and particle densities do satisfy the same 
initial conditions. To prove our point, it is enough to prove that 
the system of equations 
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(9) 



or the equivalent 



&3y 
\ An 
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(10) 



admits at least one solution. This is again trivially true, since 
this is a system of 3 equations in 4 variables. We can now as- 
sume the density An assigned and find, if the 3x3 submatrix of 
r" is invertible, the corresponding current density that satisfies 
eqn ( fTOb . If that submatrix of V is not invertible, then we will 
be able to find more than one solution for Aj, thus reinforcing 
our statement. The solution found in this way can also be made 
to fulfills the initial condition Aj(r, t = 0) = 0. 

It is also worth pointing out that the standard proof of the 
theorem for STDCDFT— , as well as that of the theorem of 
standard TD-CDFT, 29 does not require that the particle den- 
sity in the auxiliary system is equal to the particle density in 
the original one. Indeed, in the standard TD-CDFT proof one 
uses the standard continuity equation to infer the equality be- 
tween the two densities. In our proof 14 , we postulated it from 
the uniqueness of the solution of the non-linear continuity eqn 
(Q~|), with Ts(r, t) a functional of the current density. In those 
proofs one only needs to determine the n-th time derivative of 
all quantities, and only the (n + l)-th time derivative of the 
vector potential. The n-th time derivative of the particle den- 
sity is obtained from eqn (T} (or, in a closed system from the 
same equation with Tb = 0). That equation, however, does not 
contain any (n + l)-th time derivative. Therefore, the equation 
for the vector potential in the auxiliary system we use in our 
proof 14 is still a recursive relation, with a unique solution pro- 
vided the initial conditions. The theorem then guarantees the 
one-to-one correspondence between external vector potential 
and ensemble-averaged current density, leaving open the (pos- 
sibly quite difficult) task of obtaining the ensemble-averaged 
density from the current density, when the continuity equation 
is given by eqn (Q~|). 

This one-to-one mapping, implies that Tb is a functional 
of the current density j(r, t), the initial conditions and the cou- 



pling between the system and the environment. This ultimately 
implies that- by the knowledge of this functional-we could 
uniquely determine the single particle density. This is how- 
ever a common situation with DFT The theorems of DFT offer 
a solid foundation for certain calculations by guaranteeing that 
one quantity can be exactly obtained from a doppelganger sys- 
tem. The other physical quantities have to be derived from the 
first one, usually a difficult-if not hopeless-task. 

We now want to critically examine Theorem 3 in Ref.—, 
where the main result of that paper is presented: it is possi- 
ble to construct a closed non-interacting quantum system that 
mimics the dynamics of n(r, t) and j (r, t) of the real open inter- 
acting system. In the following, we will assume that Tb (r, t) 
could be written as the divergence of en(r,t)C(r,t)/m as in 
ID We want to point out that, one can in fact find many (possi- 
bly infinite) closed non-interacting quantum systems that repro- 
duce the dynamics of the exact current and particle densities. 
One can show that this is not in contrast with the general the- 
ory of TDCDFT. Indeed, we can show that the total current in 
the closed KS system is uniquely determined by the total vec- 
tor potential, given by the sum of the KS vector potential and 
C(r, t). However, the particle and the current density j(r, t) do 
not share this property, and we can find infinitely many closed 
non-interacting KS systems that reproduce these quantities. In- 
deed, in the proof of the Theorem 3 in Ref.— it is stated that 
the vector 



C(r,t) 



en(r, t) 



, / dn(r, t) ^* ->, 
d 3 r{— i-Li + V-j(r,t) ] (11) 



can be uniquely determined by eqns (Q3 and (O once a bound- 
ary condition in space that fixes an arbitrary function of time is 
assigned. 

However, the above equation is not the unique solution. For 
instance, easily fulfilling the assigned boundary condition, one 
can add the curl of an arbitrary vector, g(r, t), to the "leakage 
potential" C(r,t) and still satisfy eqn ( fTTT i. Indeed it is eas- 
ily proven that C'{r,t) = C(r,t) + (V x g(r,t))/n(r,t) still 
satisfies eqns Q] and [4] if C(r,t) satisfies the same equation. If 
we now continue with the proof and use C'(r, t) we arrive at a 
new vector potential A' KS (r, t) that gives the same n(r, t) and 
j(r, t) as the couple Af^s(r, t) and C(r, t). This ambiguity re- 
flects the fact that we are trying to mimic the effect of a scalar 
function - the term Tb(t, t) in eqnQ]- with a vector function, 
C(r, t), without imposing strict boundary conditions (BCs) on 
C(r,t). 

A simple solution to the aforementioned problem may ap- 
pear by setting V x C(r, t) = 0. However, it is important to 
realize that the imposition of certain boundary conditions on the 
dynamics of these quantities has a direct impact on the unique- 
ness of the results. For example, assuming that C(r, t) reaches 
a certain uniform limit when |r| — > oo, might be inconsistent 
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with the bath operator acting on the true many -body system. In- 
deed, certain bath operators can be strongly non-local in space, 
effectively transferring charge from one region of space to an- 
other, with the two arbitrarily far from each other. Therefore, 
the BCs on V x C (r, t) have not a clear physical origin or re- 
lation to any physical observable. Fixing their value to obtain 
one solution appears utterly arbitrary. 

Notice that a similar problem appears also in the case of the 
standard theorem of TDDFT (see for example Ref.— ) where 
the additional boundary condition n(r)V AV(r,t) when 
|r| — > oo is added to the proof. However, in this case, while in 
principle this condition is arbitrary and one may choose another 
condition, this choice is motivated by physical arguments that 
are valid for a wide range of systems. The same considerations 
instead do not apply to all the components of the vector C(r, t). 

These ambiguities derive again from the attempt to find a 
mapping between the vector potential and the particle and cur- 
rent densities, when the continuity equation is not valid. In this 
attempt, one needs to include another scalar quantity to make 
that mapping meaningful. The choice of a vector quantity like 
C (r, t) in eqn dTTb is ill-posed and leads to another ambiguity. 
On the other hand, if we are to gain any physical insight on 
the particle density, we need to obtain reliable approximations 
of the leakage potential. Instead, if one accepts that the quan- 
tity of interest is the current density, the leakage potential must 
be a functional of the current density and therefore we could 
solve the continuity equation for the particle density. For this 
reason, we can prove the following theorem (see below about 
some caveats regarding the use of a density matrix vs. a state 
vector formulation) 

Theorem A (unique mapping from open to closed system): 
Consider the dynamics of a many-body system in contact with 
an external environment. Assume the evolution of the density 
matrix is given by Eq. (O with the full many-body Hamiltonian 
([3J. Then under reasonable physical assumptions, there exists 
a closed non-interacting auxiliary many-body Kohn-Sham sys- 
tem which starting from given initial conditions for the state of 
the many-body open system, evolves according to 

id t \*Ks{t)) =H KS (t)\^ KS (t)) (12) 
with Hamiltonian 

H«,M = E( * + '£'' (f "" )' '»> 

and reproduces the dynamics of the current density of the orig- 
inal many-body open system, 

jKs{r,t)=j(r,t), (14) 

The proof of this theorem is identical to those already present 
in the literatur e 14 i 15 - 19 , we are therefore not reproducing it here 



. The only difference is that the current density of the original 
many-body open system follows the dynamics induced by the 
many -body Hamiltonian and by the external environment. 14 i 15 i ls 
The single-particle density can be obtained from eqn [T] once 
the leakage potential is written as a functional of j(r, t). The 
uniqueness of the closed quantum system then follows from 
the fact that two closed systems which share the same single- 
particle current density and initial conditions do coincide. 29 

It is also interesting to note that the possibility of studying 
the dynamics of an open system by using a TDCDFT closed 
system, has been presented in the past. For example in— the 
dynamics of a 2D electron gas coupled with the electrons con- 
fined in a ID quantum well is investigated. The quantum well, 
via the Coulomb interaction, acts as an external forcing field 
that provides energy to the 2D gas. The dynamics of the latter 
is then described by an effective vector potential which contains 
no reference to the ID well. One can show that the 2D electron 
gas will relax to a steady regime. A similar analysis has been 
performed for a ID electron gas confined in a quantum well 
by using the Vignale-Kohn functional for TDCDFT. 33 In that 
case, the system relaxes to the ground state.^ 4 ^ The KS en- 
ergy is lost to those degrees of freedom that the Vignale-Kohn 
functional does not take into account in describing the dynam- 
ics of the many-body system through some kind of common 
variable like the single-particle current density. It should not be 
surprising then that in this case the KS energy could be given 
a physical interpretation as the maximum work that could be 
extracted from the system. 35 

Finally, we want to comment on a fundamental but impor- 
tant issue. As we have discussed at length in our previous pub- 
lications,^ 4 ^ and at the beginning of this paper, one has two 
choices to describe an open quantum system: in terms of the 
density matrix or the state vector. These two formulations are 
generally equivalent, being the state vector formulation an un- 
raveling of the density matrix one. However, the density matrix 
approach suffers from two drawbacks which do not make it a 
solid starting point for a formulation of DFT for open quantum 
systems. This is due to both the possible loss of positivity of the 
density matrix if an equation of motion of such quantity is em- 
ployed with the Hamiltonian and/or bath operator(s) dependent 
on time,—, and the fact that the KS Hamiltonian does depend 
on internal degrees of freedom. Starting from the master equa- 
tion formulation of the same problem, one needs to exclude 
from the outset the possibility that the Hamiltonian of any aux- 
iliary system with different interaction potential (and hence the 
KS Hamiltonian) depends on the internal degrees of freedom. 
Otherwise, for such a system no closed density-matrix equation 
can be obtained. In other words, one needs to start from an hy- 
pothesis that constitutes part of the final thesis. It is only when 
one starts from a stochastic Schrodinger equation for the state 



t See for example the book by G. Giuliani and G. Vignale for a clear formula- 
tioni£. 
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